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Abstract 

We study the weak convergence (in the high-frequency hmit) of the frequency components 
associated with Gaussian-subordinated, spherical and isotropic random fields. In particular, 
we provide conditions for asymptotic Gaussianity and we establish a new connection with 

random walks on the hypergroup SO (3) (the dual of S'0(3)), which mirrors analogous re- 
sults previously established for fields defined on Abelian groups (see Marinucci and Peccati 
(2007)). Our work is motivated by applications to cosmological data analysis, and specifically 
by the probabilistic modelling and the statistical analysis of the Cosmic Microwave Back- 
ground radiation, which is currently at the frontier of physical research. To obtain our main 
results, we prove several fine estimates involving convolutions of the so-called Clebsch-Gordan 
coefRcients (which are elements of unitary matrices connecting reducible representations of 
5*0(3)); this allows to intepret most of our asymptotic conditions in terms of coupling of 
angular momenta in a quantum mechanical system. Part of the proofs are based on recently 
established criteria for the weak convergence of multiple Wiener-Ito integrals. 
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1 Introduction 



This paper deals with weak limit theorems involving the high-frequency components (in the 
sense of the spherical harmonics decomposition) of random fields defined on the unit sphere 
S^. Our results are motivated by a number of mathematical issues arising in connection with 
the probabilistic and statistical analysis of the Cosmic Microwave Background radiation (see 
e.g. [6]). We start by giving a description of our abstract mathematical framework, along with 
a sketch of the main results of the paper; the subsequent Section 11.21 focuses on the physical 
motivations and applications of our research. Here, and for the rest of the paper, all random 
elements are defined on a suitable probability space {Q,,J^,F). 
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1.1 General framework and outline of the main results 

We shall consider real-valued random fields {T{x) : x £ S^} enjoying the following properties: 

ET (x) = , ET2(x) < +00 and f (^x) '= f{x), (1) 

for all X € S'^ and all g G SO{3), where denotes equality in law (in the sense of stochastic 
processes). A field verifying the last relation in ([1]) is usually called isotropic or rotationally- 
invariant (in law). It is a standard result that the following spectral representation holds in the 
mean-square sense: 

00 00 / 

f{x) = j2 Ti {x) = Y.Y. ' (2) 

«=0 1=0 m=-l 

where {l^m '■ I > 0, m = —I, /} is the collection of the spherical harmonics, and the {aim} ^-re 
the associated (harmonic) Fourier coefficients. For / > 0, we also write Ci = E|a/m|^5 and we 
call the sequence {C; : / > 0} the angular power spectrum of the random field T (note that C/ 
does not depend on m - see e.g. [2]). For every I > 0, the field Ti provides the projection of T on 
the subspace of L^{§'^,dx) spanned by the class {Yi^ ■ m = —I, ...,/}. The spherical harmonics 
form an orthonormal basis of L^(§^, dx) which can be derived from the restriction to the sphere 
of harmonic polynomials. In particular, in spherical coordinates x = {9,ip) they can be written 
explicitly as: loo = 1/V'47r and 

Yim{e,v) = J^^f—^Pim{cos9)e^''^^,m>0, (3) 
y 47r (t + mj! 

Yimie, 95) = i-irYi^^ie, ^),m<0, 0<e<7r, 0<ip<2TT , (4) 
where, for / > 1 and m = 0, 1, 2, I, Pim{ ) denotes the Legendre polynomial of index /, m, i.e., 

PUx) = (-l)-(l - x^r^'^Piix) , Plix) = ^^(x^ - 1)^ (5) 

For a discussion of these and other properties of the spherical harmonics see e.g. [131 Chapter 9], 
or [281 Chapter 5] . For / > 0, the real- valued field Ti is called the Ith frequency component of T. 
The expansion ([2]) can be achieved by many different routes, for instance by a Karhunen-Loeve 
argument or by means of the stochastic Peter- Weyl theorem, see for instance [I] , [3] , [12] and [22] . 
The random harmonic coefficients {aim} appearing in ^ form a triangular array of zero-mean 
random variables, which are complex- valued for m 7^ and such that Ea^^a/'m' = (^^^ 
bar denotes complex conjugation and 5 is Kronecker's symbol; note also that aim = (— l)'"a/-m)- 
For a Gaussian random field T verifying ([U, it is trivial that the set {aim} is itself a complex- 
Gaussian array, with independent elements for m > 0. It is a simple but interesting fact that 
the converse also holds, i.e. that, under an isotropy assumption on T, the independence of the 
a/m's for m > implies Gaussianity, see [2j. Apart from this result, the behaviour of the array 
{aim} and of the projections {T;} for non-Gaussian isotropic fields is so far almost completely 
unexplored and open for research, although such objects are highly relevant for cosmological 
applications (see the next subsection). It should be stressed that the coefficients {aim} depend 
on the choice of coordinates and are not intrinsic to the field, although their law is. In this sense, 
it is sometimes physically more sound to focus on the behaviour of the sequence of projections 
{T;}, which are indeed invariant with respect to the choice of coordinates. 
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In what follows, we focus on non-Gaussian fields T that are Gaussian-subordinated, and we 
address the previous topic by studying the asymptotic behaviour of {aim} and {?]}, as I — > +00. 
Recall that T is called Gaussian-subordinated whenever T (x) = F (T (x)), where F is a suitable 
real-valued function, and T is an isotropic spherical (real) Gaussian field. In particular, our 
purpose is to establish sufficient (and sometimes, also necessary) conditions on F and on the 
law of T to have that the following two phenomena take place: (I) as / — > -|-oo, for a fixed m 
and for an appropriate sequence ti (/) (/ > \m.\), the sequence 

n (/) X aim = n (0 [ F{T (z)) Yim {z)dz, I > \m\ 

converges in law to a Gaussian random variable (real- valued for m = 0, and complex- valued for 
m 7^ 0); (II) for a suitable real- valued sequence T2 (/) (/ > 0) and for / sufficiently large, the 
finite-dimensional distributions of the field 

T2 (0 X fl (•) = T2 (/) Yj ' 

m=—l,...,l 

are close (for instance, in the sense of the Prokhorov distance - see |21] ) to those of a real 
spherical Gaussian field. Note that both results (I) and (II) can be interpreted as CLTs in the 
high-frequency (or high-resolution) sense, since they involve Gaussian approximations and are 
established by letting the frequency index I diverge to infinity. 

Our findings generalize previous results, obtained in [T7j, for fields defined on Abelian com- 
pact groups. One of our main tools is a result concerning the Gaussian approximation of multiple 
Wiener- Ito integrals established in ^21j (see also [19], [23] and [Mj)- These CLTs can be seen as 
a simplification of the combinatorial method of diagrams and cumulants (see e.g. [26] )■ These 
techniques, combined with the use of group representation theory, lead to one of the main con- 
tributions of this paper: the derivation of sufficient (or necessary and sufficient) conditions for 
(I) and (II), expressed in terms of convolutions of Clebsch-Gordan coefRcients (see e.g. \28\ 
Ch. 4]), which are the elements of unitary matrices connecting specific reducible representations 
of 50(3). Clebsch-Gordan coefficients are widely used in quantum mechanics, and admit a 
well-known interpretation in terms of probability amplitudes related to the coupling of angular 
momenta in a quantum mechanical system (see [13], [28j or Sections [3] and [6] below) . We will 
also show that many of our conditions can be alternatively restated in terms of 'bridges' of 
random walks on SO (3) (the dual of SO (3)). The definition of such random walks differs from 
the classic one given in [8j, although the two approaches can be related through the notion of 
mixed quantum state (see Section [6|). Note that an analogous connection with random walks 
on X'^ was pointed out in [T7J. 

1.2 Cosmological motivations 

The Cosmic Microwave Background radiation (hereafter CMB) can be viewed as a relic radiation 
of the Big Bang, providing maps of the primordial Universe before the formation of any of the 
current structures (approximately, 3 x 10^ years after the Big Bang); as such, it is acknowledged 
as a goldmine of information for fundamental physics. Many satellite experiments involving 
hundred of physicists throughout the world are devoted to the construction of spherical maps of 
the CMB radiation, and for pioneering work in this area G. Smoot and J. Mather were awarded 



3 



the Nobel Prize for Physics in 2006 ~ see for instance http://niap.gsfc.nasa.gov/ for more 
details. 

The crucial point is that most cosmological models imply that the CMB radiation is the 
realization of a random field {T{x) : x G S^}, verifying the three conditions in ([1]); each x E 
corresponds to a direction in which the CMB radiation is measured. The isotropic property 
can be seen as a consequence of Einstein's cosmological principle, roughly stating that, on 
sufficiently large distance scales, the Universe looks identical everywhere in space (homogeneity) 
and appears the same in every direction (isotropy). A central issue in modern cosmology relates 
therefore to the distribution of the CMB random field T, which is predicted to be (close to) 
Gaussian by some models for the dynamics at primordial epochs (for instance, by the so-called 
inEationary scenario), and non-Gaussian by other models, where fluctuations are generated by 
topological defects arising in phase transitions of a thermodynamical nature - see for instance 
[6]. Many testing procedures have been proposed to tackle this issue; in some form, they all 
rely asymptotically on the behaviour of the field at the highest frequencies (see for instance [1], 
[14j and the references therein). This is a sort of unescapable, foundational issue in Cosmology. 
By definition, the latter is a science based on a single realization, e.g. our Universe or the trace 
of its primordial structure in the form of the CMB radiation, which is observed at higher and 
higher resolutions. As such, an asymptotic theory for statistical tests is possible only in the 
sense of observations at higher and higher frequencies (smaller and smaller scales) becoming 
available as the experiments become more sophisticated. In particular, any satellite experiment 
measuring the CMB radiation can reconstruct the spherical harmonic developement appearing in 
([2]) only up to a finite frequency /max, the quantity vr/Zmax representing approximately the angular 
resolution of the experiment (the pioneering satellite COBE (1993) could reach a frequency 
Zmax ^ 20, WMAP (2003, 2006) improved this hmit to Z^ax ^ 600/800, and Planck (to be 
launched in 2008) is expected to reach l^aax — 2500/3000). In order for such procedures to yield 
consistent outcomes, one should therefore figure out what is the limiting behaviour of {T/}, 
for / >> 0, under different distributional assumptions on T. Some Monte Carlo evidence (see 
for instance [16] and the references therein) has suggested that this behaviour may be close to 
Gaussian even in circumstances where the underlying field T clearly is not. The investigation of 
this issue is necessary for rigorous inference on CMB data, and in particular for non-Gaussianity 
tests. The relevance of the asymptotic behaviour of the {Ti}, however, goes much beyond the 
issue of such tests, and relates indeed to the whole statistical analysis of CMB - which is largely 
dominated by likelihood approaches (see [7]). 

We stress by now that the results we provide cover models that are quite relevant for cosmo- 
logical applications, for instance the so called Sachs-Wolfe model, which represents the standard 
starting model for the inflationary scenario (see for instance [1], [6]). In its simplest version, 
this model implies that the CMB is a straightforward quadratic transformation of an underlying 
Gaussian fleld, i.e. 

f{x) = T{x) + fNL{T{xf-ET{xf} , xeS^ (6) 

where /tvl is a nonlinearity parameter depending on constants from particle physics and T is 
Gaussian and isotropic. As a special case, our results do allow for a complete characterization 
of the high-frequency behaviour of models such as ([6]), and in this sense they are immediately 
applicable in the cosmological literature. 
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1.3 Plan 



In Section 2 we provide some background material on isotropic random fields on the sphere. 
Section 3 is devoted to a discussion on representation theory for the group of rotations SO (3) 
and the so-called Clebsch-Gordan coefficients, which will play a crucial role in the analysis 
to follow. In Section 4 we state and prove a general CLT result for the spherical harmonics 
coefficients and the high-frequency components of a field arising from polynomial transformations 
of arbitrary order of a subordinating Gaussian process. In Section 5 we provide a more detailed 
analysis of necessary and sufficient condition for the CLT to hold in the case of quadratic 
and cubic transformations; we also highlight the connections between our conditions and the 
theory of random walks on hypergroups. The interplay with random walks on hypergroups is 
further explored in Section 6, where some comparisons with the existing literature are provided, 
and some physical interpretations of our conditions in terms of randomly interacting quantum 
particles are given. In Section 7, we turn our attention to more explicit conditions on the angular 
power spectrum, and we discuss an exponential/algebraic duality which parallels to some extent 
some earlier findings in the Abelian case. 

2 Preliminaries on Gaussian and Gaussian-subordinated isotro- 
pic fields 

As in the Introduction, we denote by the unit sphere = jxGM'^: ||x|| = l}. For every 
rotation g G SO (3) and every x G S^, the symbol gx indicates the canonical action of g on x (see 
[28\ Ch. 1], as well as Section [3] below, for further details). We will systematically write dx for 
the Lebesgue measure on §2, and we denote by (S^, dx) the class of complex- valued functions 
on §^ which are square-integrable with respect to dx. We denote by {Yi^ : I > 0, m = —I, /} 
the basis of given by spherical harmonics, as defined via ([3]) and ([4]). From now 

on, we shall denote by T = {T (x) : x G S^} a centered, real-valued and Gaussian random field 
parametrized by S^. We also suppose that T is isotropic, that is, for every g G SO (3) one has 

that T (x) T (gx), where the equality holds in the sense of finite dimensional distributions. 
To simplify the notation, we also assume that ET{xy = 1. Following e.g. [2] (but see also [3], 
[22] and [25]), one deduces from isotropy that T admits the spectral decomposition 

oo / oo 
T{x) = Yl aim;iYim (x) = T^x) , X G S^ (7) 

Z=0 m=-l 1=0 

where aim-i — /§2 T (x) Yim {x)dx (the role of the subscript "/m; 1" will be clarified in the follow- 
ing discussion), T; (x) = X]!n=-; '^im;iYim {x), and the convergence takes place in (P) for every 
fixed X, as well as in L?' (P ® dx). The next result gives a simple and very useful characterization 
of the joint law of the complex-valued array {a;^;! : / > 0, m = — /, For every z G C, the 

symbols {z) and 9 (z) indicate, respectively, the real and the imaginary part of z. 

Proposition 1 Let T he the centered, isotropic and Gaussian random field appearing in 
Then: (i) for every I > the random variable aiQ-i is real-valued, centered and Gaussian; (ii) 
for every I > 1, and every m = 1,...,/, the random variable aim-i is complex-valued and such 
thataim-i = (-!)"" ai_m;i, and moreover ¥.{^{aim-if) = E(9 (aj^-i)^) = E(afo.-^)/2 = Ci/2, for 
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some constant Ci £ [0, +00) not depending on m, and 

E(3?(az™;i)9(az^;i)) = 0; (8) 

(in) for every I > 1 and every m = —I, the random coefficient ai^-i is independent ofai'm';i 
for every I' > such that I' ^ I and every m' = By noting Cq = E(aQQ.;^), one also has 

the relation 

1=0 

The reader is referred to [2j for a proof of Proposition [1], as well as for several converse 
statements. Here, we shall only stress that formula ([9]) is a consequence of the well-known 
relation (see e.g. [28]) 

' 2/ + 1 

Yimix)Yira{y) = ^^Pi{cos{x,y)), x,y£S^ (10) 

m=—l 

where {x,y) is the angle between x and y. Observe that property ([8]) implies that 3ft (a^^;!) and 
9 {aim-i) are independent centered Gaussian random variables. Moreover, the combination of 
([8]) and point (iii) in the statement of Proposition [T] yields that ^iaim;iai'm';i) = 0, y{l,m) ^ 
{l',m'). Finally, it is also evident that points (i)-(iii) in the previous statement imply that the 
law of an isotropic Gaussian field such as T is completely characterized by its angular power 
spectrum {C/ : / > 0}. To avoid trivialities, we will always work under the following assumption: 
Assumption. The angular power spectrum {Q : ^ > 0} is such that Ci > for every /. 

Note that the results of this paper could be extended without difficulties (but at the cost of 
an heavier notation) to the case of a power spectrum such that C/ 7^ for infinitely many Ts. 
In the subsequent sections, we shall obtain high-frequency CLTs for centered isotropic spherical 
fields that are subordinated to the Gaussian field T defined above. 

Definition A {Subordinated fields). Let Lq(R, e~^'^^^dz) indicate the class of real- valued 
functions F (z) on M, which are square-integrable with respect to the measure e~^ ^^dz and 
such that J F (z) e~^ ^"^dz = 0. A (centered) random field T = {T (x) : x € S^} is said to be 
subordinated to the Gaussian field T appearing in ([2]) if there exists F € Lq{M., e~^ /"^dz) such 
that T {x) = F [F] [x), Vx S S^, where the symbol F [T] (x) stands for F (F (x)). Whenever T is 
subordinated, we will rather use the notation F [F] (x) instead of F (x), in order to emphasize 
the role of the function F. Of course, if F (z) = z, then F [F] (x) = F (x) = F (x). 

It is immediate to check that, since F is isotropic, a subordinated field F [F] (•) as in Definition 
A is necessarily isotropic. As a consequence, following again [2] or [22], one deduces that F [F] 
admits the spectral representation 

00 I 00 

F [F] {x) = Y,Y. (^) = E ^ [^]; (^) ' ^ ^ (11) 

1=0 m=-l 1=0 

with convergence in (P) (for fixed x) and in (17 x S^, P dx). Here, 

aim{F) ^ [ F[F]{y)Y;;j^dy,and (12) 

I 

F [F]i (x) ^ Yl (^) (^) • (13) 

m,=—l 
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The complex-valued array {aim (F") -1^0, m = —I, ...,/} always enjoys the following properties 

(a) -(c): (a) for every / > 0, the random variable a^o {F) is real-valued, centered and Gaussian; 

(b) for every / > 1, and every m = 1, I, the random variable aim [F) is complex-valued and 
such that 



aim{F) = {-ir ai_m{F) ; E(5R (a,„ (F)) 9 (a,„ (F))) = 
E(K (aim iF)f) = E(9 {aim iF)f) = ^aio {Ff)/2 = Q (F) /2, 



where the finite constant C; (F) > depends solely on F and /; (c) ^.{aim (F) x ai/m' (F)) = 0, 
V(/,m) 7^ {l',m'). Note that, in general, it is no longer true that ?R.{aim{F)) and Q{aim{F)) 
are independent random variables. Moreover, we state the following consequence of [21 Th. 7]: 
for every I > 1, the coefficients {aio (F) , ...,aii (F)) are stochastically independent if, and only 
if they are Gaussian. Also, E{F [T] (xf) = J^Zo ^'^i (^)- 

In the subsequent sections, a crucial role will be played by the class of Hermite polynomials. 
Recall (see e.g. |Tll p. 20]) that the sequence {Hq : q > 0} of Hermite polynomials is defined by 
the differential relation 

Hq{z) = {-iye^—e--, zeM, g > 0; (14) 

— 1 /2 

it is well-known that the sequence {{q\y'''Hq : g > 0} defines an orthonormal basis of the 
space L2(R, (27r)"^/^ e-'-"^ /"^dz). When a subordinated field has the form (for q > 2) Hq [T] (x), 
X G (that is, when F = Hq in Definition A), we will use the shorthand notation: 

(15) 
(16) 

(17) 

/ > 1, X G S^ (18) 
= -l,...,l. (19) 

To justify our notation (|15p - (|19|) . we recall that for every fixed x the random variable 
Hg [T] (x) = Hq (T(x)) is just the gth Wick power of T (x) (see for instance [H]). We con- 
clude the section with an easy Lemma, that will be used in Section [H 

Lemma 2 Let F [T] (x), x G S^, be an (isotropic) subordinated field as in Definition A. Then, 
for every I > 1 one has the following: 

1. The random field x F [T]i (x) defined in il3\} is real-valued and isotropic; 



r(«) (x) 


A_ 


Hq [T] (x) 


, X G s^ 


aim;q 


A_ 


aim {Hq) , 




r/«) (x) 


_A 


Hq [T]i (x) 


, / > 1, X G §^ 


tI'^ (-) 


A_ 


Var (t/^) 


(x) j r/^^ (x 




_A 


Cl (Hq) = 


E|azm;gp, / > 1, 



2. For every fixed x G S^, F [T]i (x) W ^^a^o (F), where the coefficient aiQ (F) is defined 



An 

according to il2\) . and consequently K{F [T]i (x)^) = ^j^Q (F); 
3. The normalized random field 



F[T]i{x) 



(21 + 1) Cl (F) 



-1/2 



F[T]i{x) (20) 
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has a covariance structure given by: for every x,y £ S^, 

E (¥W]i (^) X ¥W]i (x)) = Pi (cos {x, y)) , (21) 

where Pi (•) is the Ith Legendre polynomial defined in ^ and, as before, {x,y) is the angle 
between x and y. 

Proof. Point 1. is straightforward. To prove point 2. define (in polar coordinates) xq = 
(0,0) and use the isotropy property stated at point 1. to write 



F [T], (x) F [T], (xo) = ^irn {F) Yi^ (xo) = y ^a,o (F) , 

m=—l 

since ([3]) implies that Yim (xq) = \J {21 + 1) /iird^. Finally, to prove relation (f2T]) we use (fTOl) to 
deduce that, for every x,y £ 

2/ -I- 1 

E(F [T], (x) F [T]i (y)) = Q (F) -^Pii^os {x, y)), 

thus giving the desired conclusion (recall that P; (1) = 1). ■ 

For instance, a first consequence of Lemma [2] is that, for every q > 2, 

E(r/^) (x f) = {21 + 1) c/'^V^vr (22) 
where we used the notation introduced at ([I5|)-([l9|), so that T^''^ (x) = [{21 + 1) C'/'^V4vr]"^/^ 

The main aim of the subsequent sections is to provide an accurate solution to the following 
problems (P-I)-(P-III). 

(P-I) For a fixed q > 2, find conditions on the power spectrum {Ci : / > 0} of T, to have that 
the subordinated process 

rp{q) ^ |2"(<?) (a;) : x e §2} defined in (jlSp is such that, for every x G § , 



{2l + l)cl'^^/4Tr xT^^^x) ^ N, (23) 

/ — >+oo 

where is a centered standard Gaussian random variable. 

(P-II) Under the conditions found at (P-I), study the asymptotic behaviour, as I — > +oo, of 
the vector 



{21 + 1) Cl'^'/^TT X {Tl''> (xi) , ...,Tl'^> (xfc)j , (24) 
for every xi,...,Xfc G S^. 

(P-III) Combine (P-I) and (P-II) to study the asymptotic behaviour (in particular, the asymp- 
totic Gaussianity), as / +oo, of vectors of the type 

V(2/ + l)CKF)/4vr X {F [T\ (xi) , F [T\ (x^)) , (25) 

for every xi, x^ G and every F £ Lq (m, e^^^^'^dz 

Note that Problems (P-I)- (P-III) are stated in increasing order of generality. We observe 
also the following fact: since (PT]) holds, and since the limit of Pi ((x, y)) {I ^ +oo) does not exist 
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in general, it will not be possible to prove that the vectors in (j24p and (j25p converge in law to 
some Gaussian limit. However, by using the results developed in [21j . we will be able to establish 
conditions under which the laws of such vectors are "asymptotically close" to a sequence of k- 
dimensional Gaussian distributions. As already mentioned, to study (P-I)-(P-III) we shall 
use estimates involving the so-called Clehsch-Gordan coefficients, that are elements of unitary 
matrices connecting some reducible representations of SO (3). The definition and the analysis 
of some crucial properties of Clebsch-Gordan coefficients are the object of the next section. 



3 A primer on Clebsch-Gordan coefficients 

In this subsection, we need to review some basic representation theory results for 50(3), the 
group of rotations in M^. We refer the reader to standard textbooks (for instance, [2S| and [29]) 
for further details, as well as for any unexplained notion or definition. It should be stressed 
that most of our arguments below could be extended to general compact groups with known 
representations; however, throughout the following we shall stick to the group of rotations 50(3), 
mainly for the sake of notational simplicity. 

We recall first that each element g G 50(3) can be parametrized by the set (a,/3, 7) of 
so-called Euler angles, where 0<a<27r,0</3<7r and < 7 < 27r. In these coordinates, a 
complete set of irreducible matrix representations for 50(3) is provided by the so-called Wigner's 
D matrices D\a,l3,-f), of dimensions {21 + 1) x {21 + 1) for I = 0,1,2,... - see [28l Ch. 4] for 
an analytic expression. Here, we simply point out that the elements of D\a, /3,^) are related 
to the spherical harmonics by the relationship 



DUa, f3, 7) = (-1)™^ ^y,_„(/3, a) = ^ «) ' (26) 

from which it is not difficult to show how the usual spectral representation for random fields 
on the spheres (for instance ([2]) and ([7])) is really just the stochastic Peter- Weyl Theorem on 

= 50(3)/50(2). The reader is referred e.g. to [29] and [27] for further discussions on the 
Peter- Weyl Theorem, and to [2], [3] and [22] for several related probabilistic results. 

It follows from standard representation theory that we can exploit the family {-D'}i=o,i,,2,... to 
build alternative (reducible) representations, either by taking the tensor product family {D'i ® 
-^'^}ii,«2' o^ by considering direct sums {®\^J~\l^_i^\D''}i^^i^] these representations have dimensions 
(2/1 -|- l)(2/2 + 1) X (2^1 -|- l)(2/2 + 1) and are unitarily equivalent, whence there exists a unitary 
matrix Ci^^i^ such that 

{d^' ^D^'}= Ci^H {®\'=iU\D'} CI,,. (27) 

Here, d^i, is a {{2li + l)(2/2 + 1) x {2li + l)(2/2 + 1)} block matrix with blocks C\^(.^^^i, of 
dimensions (2/2 + 1) x {21 + 1), mi = —li, The elements of such a block are indexed by 

m2 (over rows) and m (over columns). More precisely. 



"-^/i(mi)Z2- 



{f-tlm I 



mi=-li,...,li;l=\l2-li\,...,l2+h 
m2=—l2,---,l2','m-=—l,---,l 



The Clebsch-Gordan coefRcients for 50(3) are then defined as {O/^^^^^^^l) that is, as the 
elements of the unitary matrices O/^/j (note that such matrices are real-valued, and so are the 
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^Iirnii2m2^' These coefficients were introduced in Mathematics in the XIX century, as motivated 
by the analysis of invariants in Algebraic Geometry; in the 20th century, they have gained an 
enormous importance in the quantum theory of angular momentum, where C\^_^j^^^^ represents 
the probability amplitude that two quantum particles with total angular momentum li and I2 
and momentum projections on the z-axis rrii and m2 are coupled to form a system with total 
angular momentum / and projection m (see e.g. [13]). Their use in the analysis of isotropic 
random fields is much more recent, see for instance [9] and the references therein. Explicit 
expressions for the Clebsch-Gordan coefficients of SO (3) are known, but they are in general 
hardly manageable (see e.g. [28], Section 8.2]). However, these expressions become somewhat 
neater when mi = m2 = m^ = 0, in which case one has the relations: C'^q^^q = 0, when /i+ I2+ 
I3 is odd, and, for h + I2 + I3 even, 

^1,0 ^ {-i)'-'^V2kTT[{h + I2 + h)/2]\ 

hOhO [(^^ i^y 2]l [{h -I2 + h)/2]\ [{-h +I2 + l3)mi 



ih+k-hV-ih-k+hV-i-h + h + hy- 
ih + h + h + iy. 



1/2 



The Clebsch-Gordan coefficients enjoy also a nice set of symmetry and orthogonality prop- 
erties which will play a crucial role in our results to follow (see and [T5| for an account of 
such properties). Note in particular that the Clebsch-Gordan coefficients are different from zero 
only if mi -|- m2 = m and \l2 — li\ < ^ < + ^2 (the triangle conditions). Also, from unitary 
equivalence we deduce that 

mi,m2 l,m 

Remark on Notation. Depending on the notational convenience, we write sometimes 
sums of Clebsch-Gordan coefficients without specifying the range of the indices / and/or m. 
In such cases, the range of the sums is conventionally taken to be the set of indices where 
the Clebsch-Gordan coefficients are different from zero. For instance, in (j28p one should read: 

Emi,m2 = Emi=-h,...,/i T.m2=-i2,...,i2 ^i,m = T.t^ Em=-i,...,i- Similar conventions are 
adopted (without further notice) throughout the paper. We recall also that the Clebsch-Gordan 
coefficients are equivalent, up to a normalization factor, to the Wigner's 3j coefficients, which 
are used in related works such as [T5j . 

The Clebsch-Gordan coefficients play a crucial role in the evaluation of integrals involving 
products of spherical harmonics. In particular, the so-called Gaunt integral gives 



Is- 



Ylj^mi {x)Yl2m2 {x)yim {x)dx — J (^21 + 1) — ^^'i™i'2m2^/%Z20 ' (^9) 



Relation (|29|) can be established using (|26|) . (|27|) and resorting to standard orthonormality 
properties of the elements of group representations - see [28t Expression 5.9.1.4]. More generally, 
define 

g {li,mi;...;lr,TrLr} = (x) • • • 1]^^^ (x) dx, (30) 



and call the quantity ^ {/i, mi; Z^, ""T-r} a generalized Gaunt integral. Then, iterating the 
previous argument, ior q > 3 it can be shown that (by using for instance [251 Expression 
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5.6.2.12]) 



g {li,mi; ...■,lg,mg;l,-m} (31) 



Li...Lq-2 Mi...Mg^2 




(2/i + 1)(2^2 + l)^LiO ^LiMi / 2/g + 1 



47r(2/ + 1) ^hOZ20'-^Zimii2m2Y 4^ '-^L,_20ZgO'-^L9-2M5_2Zgmq ) 

where, for g = 3, we have used the convention II^Lj^ = 0. Note that expressions such as (I3ip 
imply that the generahzed Gaunt integrals of the type (j30p are indeed real-valued. To simplify 
the expression (iST]) . let us introduce the coefficients 



Ai Ap_2 

^Xi,\2,...,\p-i;fi ^ \ \ ^\i,fj.i ^A2,At2 . . . f^^P~^'f^ 

h,m\;...;lpmp ~ / > ■" / > h,m\,l2,m2 Xi,tJ.i;l3,m3 Ap_2,Atp-2 ;Zp,mp ' 

^1=— -^1 Mp-2=— -^p-2 

These coefficients are themselves the elements of unitary matrices connecting tensor product 
and direct sum representations of SO (3), and thus it follows easily that the following orthonor- 
mality conditions hold 

Er^Ai,A2,...,Ap_i;/^~l _ ^ r „Ai,A2,...,Ap_i;/^~l _ , s 

\'^Zi,mi;...;Zpmp j ~ Z^'-' \'-"Zi,mi;...;Zpmp / " ^' W^J 

mi,...mp Ai Ap_i /x=-Ap_i 

it is important to note that due to the conditions mi + 1712 = the sums may actually vanish, 
for instance 

^^l,-^2vi^p-i;0 ^Ai,0 /-~(A2,0 ^-^p— liO /'QQ'\ 

'~^li,0;...-lpO — *-"Zi,0,Z2,0'-'Ai,0;«3,0 " ' ' '^Ap_2,0;Zp,0 • 

We have also that 

Q {li,mi; ...]lq,mq;l,-m] (34) 

E^Li,L2,---,Lq-2,l;Q ^Li,L2,...,Lq-2,l\m 
hfl;...;lqO ^li,mi;...;lqmq 

Ll...Lq-2 

Remark. The coefficients Ci^^mi' '/^m^'^ defined above admit a physical interpretation in 
terms of coupling of angular momenta in a quantum mechanical system. Consider indeed a 
system composed of p particles, say ai, ...,ap, such that ai has total angular momentum equal 

to li, and projection on the z-axis given by rrii. Then, the coefficient C^^^^mi' ■z^m^ '^ exactly 
the probability amplitude of the intersection of the following p — 1 events Ei,..., Ep_i: 
El = {«! and 02 couple to form a particle r/i with total angular momentum Ai}, E2 = {iji 
couples with to form a particle r/2 with total angular momentum A2},..., Ej = {?7j_i couples 
with Qfi+i to form a particle r]i with total angular momentum A2},..., Ep_i = {r/p-2 couples with 
ap to form a particle with total angular momentum Ap_i and projection fi on the z-axis}. 

In the sequel, we shall also need the so-called Wigner 6j (or Racah) coefficients, which are 
related to the Clebsch-Gordan by the identity (see ([281 Eq. 9.1.1.8])) 

\ h h '3 1 _ 7^/7 , N \ " f^hrrii f-ihm^ r'^z,mz, ^kme /qr\ 

I Ia k L \ '''' ^^ hmil2m2^hmil6m6'-'l3m3krn4'-'l2ni2l4m4 

m4Tne 
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where K {h, k) = [(2/3 + 1){21q + 1)] (^^lyi+h+U+h (note that the previous sum does not 
involve 1712 and ms, because of the general relation: C^*^^^^ ~ ^' whenever 7^ ^1+^2)- Although 
the Wigner's 6j coefficients play themselves a very important role in Quantum Mechanics and 
Representation Theory, for brevity's sake we avoid a full discussion on their properties; the 
interested reader can consult ([281 Ch.9]) or ([291 PP- 529-542]). 

4 High-frequency CLTs: conditions in terms of Gaunt integrals 

The aim of this section is to obtain conditions for high-frequency CLTs in terms of Gaunt 
integrals of the type ([3T]) . We start by focussing on Hermite polynomials, and then we deal with 
general subordinated fields. 

4.1 Hermite subordination 

We focus on the spherical field T^*?) {q > 2) defined in (jl5p . which is obtained by composing the 
Gaussian field T in ([2]) with the qth Hermite polynomial Hg (or, equivalently, by taking the qth. 
Wick power of the random variable T (x) for every x). Our first purpose is to characterize the 
asymptotic Gaussianity (when / +00) of the spherical harmonic coefficients {aim-q} defined 
in (USD. 



Theorem 3 Fix q > 2. 

1. For every I > 1, the positive constant C^'^'^ in lilS^) (which does not depend on m) equals 
the quantity 

^! ^ ••• ^ QiQa ■ ■ - Cig \G{li,mi; ...■,lg,mg;l, -m}f (36) 

h,mi lq,mq 

h,...,lq = U=l ) Ll...Lq-2 

for every m = where the (generalized) Gaunt integral G {•} is defined via [3U\) . 

2. Fix m 7^ 0. As I ^ +00, the following two conditions (A) and (B) are equivalent: (A) 

(ci'^Y'^'xai^./^N + tN', (38) 

where N, N' ^ M (0, 1/2) are independent; (B) for every p = + 1, q — 1, if q — 1 is even, 
and every p = q/2, ...,q — 1 if q — 1 is odd 



"Ijl "2(q-p) j2(9-p) 



li.nii Ip.mp 



xg{li,mi;...;lp,mp;ji,ni;...;jg_p,ng^p;l,-m} x (39) 



xQ{li,mi; /p, m^; jg_p+i, n^.p+i; J2(q-p)> ^2(g-p); -1^} 







3. Let N be a centered Gaussian random variable with unitary variance. As I — > +cxd, the 
CLT 

(Cl'^Yy^xaio^q'^N (40) 
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takes place if, and only if, the asymptotic condition i39\) holds for m = and for every p = 
+ 1, q — I, if q — 1 is even, and every p = q/2, q — I if q — 1 is odd. 

Proof. Consider a standard Brownian motion W = {Wt : t G [0, 1]}, and denote by L^. ([0, 1]) 
= ([0, 1] , dX) the class of complex- valued and square integrable functions on [0, 1], with re- 
spect to the Lebesgue measure d\. Now select a complex-valued family {gim '■ I > 0, — I 1^ rn < 1} CI 

([0, 1]) with the following five properties: (1) gio is real for every I > 0, (2) gim. = f— 1 
(3) / gimW 



9lm — [~^) 9l~rm 

dX = 0, V (/, m) ^ {I', m'), (4) / ^ (gim) 9 {gim) dX = 0, (5) / 3f? {gimf dA = / 9 {girrf dX 



= J gf^dX/l = Ci/2, where {Q : I > 0} is the power spectrum of the Gaussian field T. According 
to Proposition [H the following identity in law holds: 

: / > 0, -l<m<l}^= {h (gim) ■.l>0, -l<m<l}, 

where h{gim) = Jq 9imdW = 'St{gim)dW + i '^{gim)dW is the usual (complex-valued) 
Wiener-Ito integral of gim with respect to W. From this last relation, it also follows that, in 
the sense of stochastic processes, T (x) h (YI'i^o Ylm=-i 9imYim {x)^ (note that the function 
z 1-^ m9im{z)Yi„i (x) is real-valucd for every fixed x e and with norm equal to 1). Now 
define ^^^^([0,1]*^) to be the class of complex-valued and symmetric functions on [0,1]'', that 
are square-integrable with respect to Lebesgue measure. For every / G ^ ([^' ■'-]'')' define 
Iq (/) = Iq {^[f))+\Iq (9(/)) to be the multiple Wiener-Ito integral, of order g, of / with respect 
to the Brownian motion W (see e.g. |181 Ch. 1], or [H]). From the previous discussion it follows 
that, for every q>2, 



r(«) {x)=Hq{T{x))^= I 



9lmyim{^ 



. /=0 m=-l 



(41) 



where the equality in law holds in the sense of finite dimensional distributions and, for every 



/ G ([0, 1]), we use the notation f®"^ (ai, ...,aq) = / (ai) 



X / (aq) . Note that, to obtain 



the last equality in (j4ip . we used the well-known relation (see e.g. [Uj): for every real- valued 



/ G L|([0, 1]) such that 



i|([o,i]) 



1, it holds that Hq [h (/)] 



Now set h 



(-l)"'Eh,mi • • • Ez„m,5/imi " " " 9l,ni, Gih^mi; lq,mq;l,-m}, so that 



ai 



law 



52 



1=0 m=-l 



Yim {x)dx = Iq 



(g) 

Lm 



(42) 



so that (I36p follows immediately from the well-known isometry relation: 



E 



''Lm 



L2([0,1P) 



(to obtain (j42p we interchanged stochastic and deterministic integration, by means of a standard 
stochastic Fubini argument). To prove that (f37]) is equal to ((36]) . observe first that (f32|) yields 
that 



E\ ^ ^Li,L2,.--, 
/ J li,m\;...;l 



h ,1^1 



5r\..5r 



mi- 
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(the RHS of the previous expression does not depend on m). Then, use to deduce that 



^ ••• ^ g {li,mi;...;lg,mq;l,-my 



mi =— /i 



An 



21 + 



in 



2L + 1 



.1=1 



47r 



Li...Lq_2 



Ll,L2,...,Lq-2 



,i;oy 



This proves Point 1 in the statement. To prove Point 2, recall that, according to [n\ Proposition 
6], relation (|38p holds if, and only if. 



(<?)^-2 



0, 



L2([0,l]2('?-P)) 

for every p = l,...,q — 1, where the complex- valued (and not necessarily symmetric) function 



^/ m ®p ^im (which is an element of L'^{[0, 1]'^*-'^ ^-')) is defined as the contraction 



^'fX®P 4m («l'-'«2(g-p)) 
^[0,11'' ' 



(43) 



for every (oi, a2(g_p)) G [0,1]^*-'' where dxp is the Lebesgue measure on [0,1]^. Since, 



trivially, \\h!f^ ^/ mll^ ~ ll^i m ^q-p ^1 mll^ stress that, in the last equality, the first norm 
is taken in -L^([0, l]^'-'^"^'*), whereas the second is in L^([0, 1]^^) ), one deduces that it is sufficient 

to check that the norm of 4m®p^i m asymptotically negligeable for every p = + 1, q — 1, 
if g — 1 is even, and every p = q/2, ...,q — 1 if g — 1 is odd. It follows that the result is proved 



M 112 



,(<?) 112 



once it is shown that, for every p in such range, the norm \\h 



{<?) 



Lm "'Lm 



(9) ||2 



equals the multiple 

sum appearing in ([39|) . To see this, use (|13|) to deduce that (recall that Gaunt integrals are 
real- valued) 



41 ®P 4,1 (ai>->«2{g-p)) 

~ ■ ■ ■ djmi ■ ■ ■ gjq-pnq-pdjq-p+inq-p+l ' ' ' 9j2(g-p) "•2(9-p) 

^ • • • ^ Gi • • • ...■,lp,mp;ji,ni; jg^p,nq^p;l, -m} 

/l,mi lp,mp 

G {li,mi; lp,mp; jq^p^i,nq^p^i; ■■■]32(q~p),n2(q~p)]l, -m] , 

and the result is obtained by using the orthogonality properties of the gjnS. Point 3 in the 
statement is proved in exactly the same way, by first observing that a/o;g is a real-valued random 
variable, and then by applying Theorem 1 in [19]. ■ 



Remark. One has the relation E 



r(9) (a 



q\ [E{T{xf]]''. This equality can be proved 



in two ways: (i) by exploiting the representation of T*^*^-* ( multiple Wiener-Ito integral, 
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Yli T^^z*'''^ ^'^d the by expanding C^'^' according 



or (ii) by using the equahty E T^'^^ (x) 
to Theorem [3l so that one can apply the orthogonahty relations ([32 



Now recall that, according to part 2 of Lemma El T^"^^ [Xf — 
(j22p holds. This gives immediately a first (exhaustive) solution to Problem (P-I), as stated in 
Section [2l 

Corollary 4 For every q >2 the following conditions are equivalent: 



law 



^^^aiQ-q, SO that relation 



1. The CLT (2^ holds for every x E . 

The asymptotic relation i39\) takes pla 
q — 1 is even, and every p = q/2, q — 1 if q — 1 is odd. 



2. The asymptotic relation /i39\) takes place for m = and for every p = ^^J- + 1, ...,q — 1, if 



To deal with Problem (P-II) of Section [21 we recall the notation T^f' (indicating the lih. 
normalized frequency component of T^'i^) introduced in ([ISp . We also introduce (for every / > 1) 
the normalized Ith frequency component of the Gaussian field T, which is defined as 



Ti (x) = 



Tijx) 
Var{Ti (a;))V2 



Ti{x) 



(44) 



According to Lemma[2](in the special case F {z) = z), Ti is a real-valued, isotropic, centered 
and Gaussian field. Moreover, one has that E[Ti {x)Ti (y)] = E[T[''^ (x)T[''^ (y)] = Pi ((x,y)), 
for every q >2 and every Z > 1. The next result - which gives an exhaustive solution to Problem 
(P-II) - states that, whenever Condition 1 (or, equivalently, Condition 2) in the statement of 
Corollary [H is verified (and without any additional assumption), the "distance" between the 

finite dimensional distributions of the normalized field r[''^ and those of Ti converge to zero. 
For every k > 1, we denote by P(M'^) the class of all probability measures on M.^. We say that a 
metric 7 (•, •) metrizes the weak convergence on P{Mf') whenever the following double implication 
holds for every Q e P(R'=) and every {Q; : / > 1} C P{R'') (as / +oo): 7 (Q/, Q) ^ if, and 
only if, Qi converges weakly to Q. The quantity ^{P,Q) is sometimes called the '-^-distance 
between P and Q. 

Theorem 5 Let q>2 be fixed, and suppose that Condition 1 (or 2) of Corollary^is satisfied. 
1. For every k>l, every xi, ...,Xf^ G and every compact subset M C M*^, 



sup 

(Ai,...,Afc)eM 



E 



eiE-=iA,T(«'(x,) 



E 



(45) 



2. Fixxi,...,Xk and denote byC\Ti ;xi,...,Xfcj and £ (Ti;xi, ...,Xk) (I > respectively, 

the law of ^T^'^'* (xi) , ...,t|''^ {xk)^ and the law of (Ti (xi) , ...,Ti (x^)). For every metric 
7 (•, •) on P(M^) such that 7 (•, •) metrizes the weak convergence, it holds that 



lim 7 (£ (r;^^;xi, ...,Xk 

\-oo 



,C (Ti;xi,...,Xk 
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Proof. The crucial point is that the spherical field x T^f^ {x) lives in the qi\i Wiener chaos 
associated with the Gaussian space generated by T. By using this fact, and by arguing as in 
the proof of Theorem [3l one can show that the vector {T^^ (xi) , t|'^^ {xk)) is indeed equal in 
law to a vector of multiple Wiener-Ito integrals, of order g, with respect to a Brownian motion. 
Since each element of this vector converges in law to a standard Gaussian random variable, one 
can directly apply Theorem 1 and Proposition 2 in [21j to achieve the desired conclusion (see 
also [21, Proposition 5]). ■ 

4.2 General subordination 

We now give a solution to Problem (P-III), as stated at the end of Section [21 where F is a 
general real-valued function belonging to the class Lq (^,e~^^ ^'^dx^ . The function F admits a 
unique representation of the form 

F{z) = Y,''-^H,{z), zGR, (46) 

where the Hermite polynomials Hq are given by (|14p and the real coefficients Cq (F), q = 1, 2..., 
are such that 

< +00 . (47) 

As a consequence, for every / > 0, the frequency component F [T]^ (x) defined in (jl3p admits 
the expansion 

F[n{x) = f2'-^Ti'^\x), xE§^ (48) 

9=1 ^' 

where the series converges in (P) for every fixed x. Formula (08]) combined with Lemma [2] 
yields also that 



E{F [T]i (x) F [T], (y)) = (cos (x, y)) c/^^, 

where C^'^^ is given by p9]) or, equivalently, by ([37|) . The next result characterizes the asymptotic 



Gaussianity of i^'-subordinated spherical random fields. Recall the definition of F [T]^ given in 
pop . The proof is standard, and therefore omitted (it can be obtained e.g. along the lines of 
[ini Th. 4]). 



Theorem 6 Suppose that the following relations hold 

21+1 



1. For every q > 1, limi^+oo ^ [^) Cl"' /E{F [T], (x)^) ^ a'q G [0, +oo); 
2- En.>i{cAF)/Ql}'^',=crHF) <+^; 

3. For every q > 2, the asymptotic relation l^39\) takes place for m = and for every p 
+ 1, q — 1, if q — I is even, and every p = q/2, q — 1 if q — 1 is odd; 
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4. lim,^+oo lim, {21 + 1) ZT=p+i {^ycl'^ = 0- 

Then, for every k >1, every xi, ...,Xk € and every compact M C 

i<T2(i^)'/'E'=iA,r,(x,)" 



sup 

(Ai,...,Afc)eM 



E 



eiE-=iA,-F[T];(x,) 



-E 



where we used the notation ^44\ )- particular, the last asymptotic relation implies that, for 
every 7(-, •) metrizing the weak convergence on P(M'^), the ^-distance between 



(F[r],(xi),...,F[r],(xfc)) 

and (T^ {F)^^"^ [Ti {xi) , T; (x^)) converges to zero as I +oo. 

Remark. A sufficient condition, ensuring tliat points 1 and 3 in the statement of Tlieorem 
[6] are verified, is tlie following: there exist constants p{q) > such that (a) (2/ + 1) C^'^^ < p (q) 

for every g > 1 and every /, and (b) Yl'^i (j^^^^ p{^) < +oo. 

5 Explicit sufficient conditions: convolutions and random walks 

In this section, we further explicit the conditions for the CLTs proved in Section S] for the 
(Hermite) frequency components T^'^'^ , I >0- In particular, we shall establish sufficient conditions 
that are more directly linked to primitive assumptions on the behaviour of the angular power 
spectrum {C; : Z > 0}. The results of Section [Ql and Section [5131 cover, respectively, the case 
q = 2 and q = 3. Section 15.41 contains some partial findings for the case of a general q, as well 
as several conjectures. These results will be used in Section [7] to deduce explicit conditions on 
the rate of decay of the angular power spectrum {Ci : I > 0}. 

Our analysis is inspired by the following result, which is a particular case of the statements 
contained in [17, Section 3], concerning fields on Abelian groups. Consider indeed a centered 
real-valued Gaussian field V = {V (9) : 9 €z T} defined on the torus T = [0, 2it) (that we regard 
as an Abelian compact group with group operation given by xy = (x + y) mod(27r)). We 

suppose that the law of V is isotropic, i.e. that V (9) ^= V {x9) (in the sense of stochastic 
processes) for every x G T, and also EV [9)^ = 1. We denote by V (9) = X^^g^ '^^/^''^ the Fourier 
decomposition of V, and we write = E|a;|^ (note that = F^^). Fix q > 2, and consider 
the Hermite-subordinated field Hq [V] (9) = Hg {V (9)), where q is the qth Hermite polynomial. 
The Fourier decomposition of Hg [V] is Hg [V] (9) = ^i^zzof'' e^^^ . We write N,N' to indicate 
a pair of independent centered Gaussian random variables with common variance equal to 1/2: 
in [T7] it is proved that to have the high-frequency CLT 

a|^^ LHg[V](9)e-'^^d9 law , 

it is necessary and sufRcient that, for every p = 1, ...,q — 1, 

lim supP [[/p =j\Ug = l]=0, (50) 
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where {Un n > 0} is the random walk on Z whose law is given by [/q = 

F[Un+l=j I Un = k]=Tj_^. 

Note that the law of the random variable Un has trivially the form of a convolution of the coef- 
ficients (see also the discussion below). The correspondence between (j49p and the "random 
walk bridge" ([50|) has been used in [T7j to establish explicit conditions on the power spectrum 
{r^"} to have that (09]) holds. 

In what follows, we shall unveil (and apply) an analogous connection between the CLTs 
proved in Section U] and some specific convolutions and random walks on SO (3). 



5.1 Convolutions on SO (3) 

In the light of Part 3 of Theorem[3]and by CorollaryHl we will focus on the sequence {a/o;g : ^ > 0} 
(see (fT6|) ). whose behaviour as I +oo yields an asymptotic characterization of the fields T^^'^^ (•) 
defined in (jl7p . A crucial point is the simple fact that the numerator of (j39p . for m = 0, can 
be developed as a multiple sum involving products of four generalized Gaunt integrals, so that, 
by ([5T|) , the asymptotic expressions appearing in Theorem [3] can be studied by means of the 
properties of linear combinations of products of Clebsch-Gordan coefficients. As anticipated, a 
very efficient tool for our analysis will be the use of convolutions on N, that we endow with an 
hypergroup structure isomorphic to SO (3), i.e. the dual of SO (3). This will be the object of 
the subsequent discussion. 

Prom now on, and for the rest of the section, we shall fix a sequence {Ci : Z > 0}, representing 
the angular power spectrum of an isotropic centered, normalized Gaussian field T over S^, as in 
Section [2j Whenever convenient we shall write 

^ {21 + l)Ci , l>0, (51) 

so that, for I > 1 and up to the constant l/47r, the parameter Ti represents the variance of 
the projection of the Gaussian field T in ([2]) on the frequency /: indeed, according to Lemma 
[21 Var(Ti) = Ti/4tt. Also, we define the following convolutions of the coefficients Ti (in the 
following expressions, the sums over indices k, Li ... range implicitly from to +cxd): 

r2,i=Y.ri,Ti,{cl%^o? , (52) 

ll,l2 

^3,i=^^2M'^hiCLiOhof= ^ Ti^Ti,Ji.^^{C['^l\ff^i^Q)'^, ... (53) 

f,,= £r,_i,,^_,r,,(ciV,oM^'= ^r,,.r,,''^ {cl^Xfy^'y. (54) 

Ll,lq l\...lq Ll..Xq_2 

We stress that the equalities in formulae ([53]) and ([M]) are consequences of ([33]) . It will be also 
convenient to define a *-convolution of order p > 2 as: 



h Ip L\...Lp^2 

I2 Ip Li...Lp-2 
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Note that the number of sums following the equalities in formula (j55p is p — 1: however, we 
choose to keep the symbol p to denote *-convolutions, since it is consistent with the probabilistic 
representations given in formulae ()59p and (j60p below. The above *-convolution has the following 
property: for every p = 2, ...,q 

^fg+i-pj,f*pi.i^ = fg^i , and, in particular, ^ThT*i.i^ = fg^i . 
h h 

The *-convolution of order 2 can be written more explicitly as 

r2,/;Z,=Er^^(^'iVo)'- (56) 
h 

Remarks. (1) [Probabilistic interpretation of the convolutions) Write first = X^z 
(plainly, in our framework = 47r, but the following discussion applies to coefficients {F/} 
such that F* > is arbitrary) so that / i — > F//F=k defines a probability on N. The second 
orthonormality relation in (j28p implies that, for fixed li,l2, the application / i — > (^i%/2o)^ ^ 
probability on N. Now define the law of a (homogeneous) Markov chain {Z„ : n > 1} as follows: 



lo 

It is clear that P{^<j = 1} = ^q,l/ (r*)^, and also, for p >2, 

f 



P{Zi = = r//F, (57) 
F, 



(F* 



P,i.h _ ]p{Zp = l\ Zi = li} (59) 



^p,uS'^+i-P,h ^ r {{Zg = I) n {Zg^.^p = h)} {q>p-i). (60) 

The following quantity will be crucial in the subsequent sections: 

r-^i-p,,.Arp,A ^ r;^i-p,^;Arp.A ^^^^^^^ , ^ ^^^^y^ (g^^ 

observe that the last relation in (I6ip derives from 

r*+i-p,i;x/ (r.y-'' = n (Zg+i-p = I) I (Zi = A)} = ¥{{Zg = A) I {Zp = /)} , 

where the last equality is a consequence of the homogeneity of Z. Note also that we can identify 
each natural number I > with an irreducible representation of SO (3). It follows that the formal 
addition h + l2 = EiKCl%i^ 

q)^ may be used to endow SO (3) with an hypergroup structure. 
In this sense, we can interpret the chain {Z„ : n > 1} as a random walk on the hypergroup 
SO (3), in a spirit similar to [8]. In Section [6l we will discuss a physical interpretation of these 
convolutions and establish a precise connection between the objects introduced in this section 
and the notion of convolution appearing in [8j. 
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(2) {A comparison with the Abelian case) In [17], where we dealt with similar problems in 
the case of homogenous spaces of Abelian groups, we used extensively convolutions over Z. This 
kind of convolutions, that we note Arq,i (o' > 2, Z G Z) are obtained as in (f52]) - (f56]) . by taking 
sums over Z (instead than over N) and by replacing the Clebsch-Gordan symbols {C^^qi q)^ 
with the indicator 11-^^1^=1. Note that these indicator functions do indeed provide the Clebsch- 
Gordan coefficients associated with the irreducible representations of the 1-dimensional torus 
T = [0, 27r), regarded as a compact Abelian group with group operation xy = {x + y) (mod(27r)) 
(this is equivalent to the trivial relation gi'i^e''^^ = Ylii '^h+h=i^^^^ = gi('i+'2)a::)_ Note also that in 
the Abelian case one has A^*p = A^p,i-h- Also, if Ti = Vj , where {TJ} is the power spectrum 

of the Gaussian field F on T appearing in (jl9]) . one has that a^^i = \Uq = I], where {Un} is 
the random walk given in (jSOp . 



5.2 The case q = 2 

In this subsection, we provide a sufficient condition on the spectrum {Ci : / > 0} (or, equiva- 
lently, on {F/ : / > 0}, as defined in ([5T|) ) to have the CLT ([lO]) in the quadratic case q = 2. 
This condition is stated in Proposition [51 and is obtained via some preliminary (technical) 
computations and lemmas. 

According to Part 3 of Theorem O to deal with (j4U|) we shall find sufficient conditions to 
have that ([39]) takes place for m = 0, g = 2 and p = 1. From ([37[) we deduce 



?^(2) _ 9 J (2/i + 1)(2^2 + 1) ^ ^ .wo ^2 

2^ 47r(2Z + 1) ^'i^'2((-ii0/20. 

«1,«2=0 ^ ' 



(62) 



On the other hand, the multiple sums appearing in the numerator of (|39p become ((7 = 2, p = 1) 



jl,ni J2,n2 



^ Ci^Q{li,mi-ji,ni-l, -m}g{li,mi;j2,n2;l,-m} 



h,mi 



I V F F 



^ h'^hmijini'^h0ji0'^limij2n2'^hQj20 

I, mi 

E-p /^Im film, f-ilO 

ll'^limijini'^llQjl0'-^hmij2n2'-'li0j20 



[471(21 + 1)Y . ^ 

1 \" \ " T T T r r''o r'W r^io ^ 

\4'k{21 + 2^2^^^^^ '1 h'-'hOno'-'hOj20^i20ji0^i 

) \ ^ film film film filr. 

\ Z-^ ^Iimijini^limij2n2^l2m2jini^l2 



m 

l20j20 



Im 
m2i2«2 



(63) 



^ nin2mim2 
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Now, from [28l Eq. 8.7.4.20] we deduce that 



E\ ^ /~ilm /~ilm f-^lm 

(_l)/3^(2, + l)(2/ + l)(C/-, 
^(2s + l)(2; + l)(Q', 



im \2 



h jl M I ^1 J2 I 

I S I2 ] \ I S I2 



where P = h + ji + h + h-, and we used the Wigner Qj symbols, as defined in (j35|) . The last 
equality follows because the quantity h+ji +12 + 32 +2/ must be necessarily even, and therefore 
f3 must be even as well. It should be noted that the role of the pairs and (/i,mi) is 

perfectly symmetric, so we obtain also 

E\ ^ ^Im /~ilm f^lm /~ilm 



nin2 m\m2 



^(2s + l)(2/ + l)(C, 



Im \2 
ImsO) 



I S j2 \ \ I S J2 



whence 



X:(2. + l)(2i + l)(C;™,)^{ 'I 'l I } (64) 



Lemma 7 For all integers l,li,l2, ji, j2 it holds that, for some positive constant c, 
J](2s + l)(2/ + l)(Q'o%)' 

s 

< cmax 



^1 Jl M / ^1 •?2 I 
I S I2 \ \ I S I2 



1 1 1 



^/W+i ¥W+i' ^/W+^ ^/W+i 

Proof. Assume without loss of generality ji,j2 > otherwise we focus on (I65p rather than 
61) . For a G (0, 1), we have that 



^(25 + i)(2/ + i)(q 



10 \2 f h jl M / h J2 I 



lOsO) 



I 5/2(1^ si 



1-2 



s>l'l ^ ^ ■> ^ 
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- ' ^ '(2i + l)(2ii + l) fJ^^/TJi ^' ' ' 

where the last equahty has been obtained by setting a = 2/5. The second last step follows 
because ji,i2 > h^h implies ji, j2 > in view of the triangle inequalities li + j2 > ^; 

also, we used the inequality {maXs>i!j«(Cjg^g)^} < see Lemma [8] below. The bound with 

I2 can be obtained by exploiting the symmetries of the 6j coefficients; in particular, we recall 
that (see ([Mi Eq. 9.4.2.2])) 

h ji I ] ^ i I ji ^2 1 f h ji I 
I s I2 ) \ h s I ) \ I s li 

Remark. The bound provided in Lemma ([7|) is sufficient for our purposes below and we 
did not investigate its efficiency in detail. We remark, however, by setting ji = j2 = 0, we have 
explicitly (see [Ml Eq. 8.5.1.2]) 

E/^lm f-ilm /-lira /-ilm 

nin2mim2 

E/^lni /~tlm /^Im /-ilm -i 
'-"«imi00*-'/imi00'-"«2m200'-"«2m200 — ^ ' 

mim2 



Lemma S As h ^ +00, C/^^^q = 0{^) 



Proof. Unless the triangle condition 21 > li is satisfied, the Clebsch-Gordan coefficient is 
identically zero and the bound is trivial. Now recall that 

io V2lTT[(2/ + /i)/2]! pi!(2/-/i)!^i! )^/^ 

[^i/2]![(2Z-/i)/2]![/i/2]!l (2/ + /1 + I)! J ' 

For sequences {a/} and {bi}, write a; hi when both o/ = 0{hi) and hi = 0{ai) hold true. From 
Stirling's formula 



,0 ^/2m [(2/ + l^)/2fl+h)/2+l/2 j- ^2^,+!^^^ _ ^^^(2^-/0+1/2 

[/l/2]'^+l [(2/ -h + l)/2](2'-'i)/2+l/2 j (2/ + + 1)2^1+3/2 

_ V2lTl(2/ + /i)(^'+'^)/^+V^ f /f^+^(2/-/i)(2'-'^)+V2 |'^' 

~Zll+l(2/-/i + l)(2'-'i)/2+l/2 I (21 + h + j 

V^TTT 1 ^/2rTT _ , 1 , 



1/2 



;V2(2/_ /i + i)i/4 (2i+li + l)V4 - ;i/2(2/_/, + i)i/4 
■ 

We can finally state a sufficient condition for the CLT (j4U|) in the case g = 2. 
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Proposition 9 For q = 2, a sufficient condition for the CLT is the following asymptotic 
relation ^ 

lim sup ^';^'^^^;'^°\^ = lim supP{Zi = /i|Z2=/2}=0, (66) 

where the {F/} are given by 151\) and {Zi} is the Markov chain defined in formulae (51) and 
Ml- 

Proof. In the sequel, we shall use repeatedly the trivial inequality 



E 

i=o 



< 



j=0 



j=0 



(67) 



j:aj<bj j-aj>bj 

which holds for arbitrary n and real vectors {aj}, {6j} and {cj}. In view of LemmaEl by using a 
generalized Cauchy-Schwartz inequality, (j67p and symmetry considerations, we obtain that the 
expression (|63l) is such that 



< 



< 



< 



[47r(2/ + 1)]^ 



[47r(2/ + 1)]2 ^ 



r<io r^io wo riio 

'-^«l0ii0'-^«l0i20*-^i20il0^«20j20 



^io r<^o 

<-^«l0ii0'-'/l0j20'-'«20jl0'-^i20i20 



'2ii 

The last expression is less than 
8[7r(2/ + 1)]2 



8[7r(2/ + 1)]2 ^ {'^'iVio} ^^'^^^'^ {^'iVao} 

^wX^^'a^Jl {^;2Vo} X]^'2-'"j2 {'^«2V-2o} ■ 



Z]^7# {'^'iVio} Z]^'i^i2 {'^/iVo} • 

hji V Jl iij2 



Now 



J max 

h<3* 



«i>0 



;i,ii>o 



It follows that 



< 



^hji 5/-2 X'^hOjiOj 



< 2 



J* maxjj<j* 



Z]^',Z2=o ^^i-'^^2(C'/i%2o)^ 
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+ 



(68) 



Now fix e > 0. Under (j66p we have that, for any fixed and positive number > 1/e, 



lim 



J maxji<j* 



+ 



< f hm sup ■ 



Because e is arbitrary, the proof is concluded. ■ 

Remark. Note that, using (I54p and ()56p . condition ()66|) becomes 



nm sup 



. 



(69) 



Note also that if, in the convolutions (jMj) . one replaces each squared Clebsch-Gordan coefficient 
(C'J^q^^q)^ by the indicator l;^_|_/2=z and extends the sums over Z, one obtains the relation 



hm sup 



0. 



(70) 



In particular, when {F;} = {rj^} (the power spectrum of the field y on T given in (I49p ) it is 
not difficult to show that formula (j70p gives exactly the asymptotic (necessary and sufficient) 
condition (f50]) . 

5.3 The case g = 3 

Our results for q = 3 closely mirrors the conditions we derived in the previous subsection. 
Proposition 10 A sufficient condition for the CLT when q = 3 is 



lim sup 



lim sup 



l—>oo 



0, and 



. 



(71) 



(72) 



Remark. In the light of (j54p -(j56p and of the definition of the random walk Z given in ([5 
and (j58p . it is not difficult to see that (j7ip can be rewritten as 



hm sup — = lim sup — 

= lim supP [Z2 = A I Z3 = /] = 0. 

1^00 X 



(73) 
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Likewise, one obtains that (j72p is equivalent to 



iim sup r 

Z— >oo ji T T T ) y^LilQ \' 

lim supP[Zi = ji I Z3 = /] = . 



(74) 



It should be noted that the two conditions (j73p and (j74p can be written compactly as 

. (75) 



r • r* 



Iim max sup 

i^oo 9=1,2 „ „ ( ^LilO 1" 

Relation (I75p parallels once again analogous conditions established for stationary fields on a 
torus - see [17J. 

Proof of Proposition [TOl In view of Part 3 of Theorem [3l we shall focus on the asymptotic 
negligeability of the ratio appearing in ([39]) . in the case where q = 2> and p = 2. As before, the 
denominator of ()39p is proportional to 



1 1=1 ) Li 



(76) 



^1,^2, '3 
1 



(2/ + 1 



^'l,'2,^3 Li 

On the other hand, the numerator is proportional to 



(2/ + 



]1,J2 "1,"2 



L2 



^h^h ^hQl2QjiO'^hmil2m2hni ^h0l20j20^hmil2m2j2n2 
Ii,l2,ini,m2 Li 



{21 + 1)2 



J1,J2 "1,«2 



^ '1^ '2 "^«i0i20ii0 "-'iimii2m2^l^^iil"i 
Ii,l2,mi,m2 Li Ml 

E^L2lO Sr^ (^L2M2 (^Im 
^/l0«20i20 '-'«imi«2m2^i2M2j2n2 
L2 M2 



(77) 
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This last expression equals in turn 



{21 + 1)2 ^ 



(2r+i) 



LiMijini 
Ml L2 



M2 



4_ n2 X] X] ^Ji^. 

Jl,j2 ™li"21 



^ '1-^ '2 '-'«i0/20ii0'-'«i0«20j20 LiMijim'^ LiMij2n2 



h,l2=0 



Li 



Ml 



{21 + 1) 



jl,i2 "l.'^21 



M2 



EP P V^^LiZO V^^im V^^L2«0 V^Wm i:L2):- 

ii^ h / ^'-'hOhOhO/ LiWlihni / J^hQUQ^o^) / >'-^I/2M2.72n2"Li "Mi 

/l,«2 -J-i Ml L2 M2 



and we can use the same argument as for q = 2. More precisely, one can write 

2 

iLilQ \ ^ /^(m 
ii,«2=l ii Afi 



E ''^'i-'^'a E ^iiW20ji0^«i0i20i20 '-^LiMiiini'-'LiAfij2n2 



Ep p ^Li«0 piLilO ^L2l0 (^L2l0 

^h-'-^h 2^ '-"ii0«20ji0'-"«i0i20i20'^«30U0ii0'-'«30«4i 



«30«40j20 



;i.../4 



L1L2 



E/^im (-lira (-ilra (-~<lra 

LiMijim'^ LiMij2n2'-' L2M2jini^ L2M2j2n2 



M1M2 



EP p ^Li«0 ^Li«0 (-^L2l0 (^L2l0 

^Ii---^Ia '-"iiOZ20iio'-"iiOi20i2o'^«30UOiiO*-"Z30«40j20 

Ii...Ia L1L2 

-,10 \2 { Li ji I \{ Li j2 I 



(-i)^X(2.+i)(2/+i)(q 



lOsa) 



I s -L2 J I Z s L2 



where C = -^1 + Ji + -^2 + j2) and (jTHj) equals 



EP p ^LiZO ^Li«0 ^L2«0 f^L2lO 

i-h-'-i-U '-^«i0«20ji0'-'ii0«20j20'-'i30U0ii0'-^Z30U0i20 



il..J4 L1L2 

{-lfY}2s+l){2l + l){ClZo) 

s 

From d?]) we now obtain that the last expression is bounded by 



/ s L2 \ \ I s L2 



EP p \" ^LilO (yLilO (-1L2IO 
i«i---Lu 2^ "-^Zio;20iio'-^/iO/20j2o'-^i30Z40iio'-"/ 



L2«0 



;i..j4 



L1L2 



+ E T/.-ru E 



LiiO ^LiiO ^L2«0 ^L2i0 



1 



«l0«20jl0^«l0«20j20^«30«40jl0^«30U0i20 5 



;i..j4 



L1L2 



Jl 
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whence all the terms are bounded by 



3132 h-.-h L1L2 



1 



_ rilO ^L20 ^/O ^-L20 

X *-^Ll0i20'-^«30«40'^i20il0'-^Z30Z40'-"i20i20"i7^ 



3132 h—U L1L2 



X C 



10 



riL20 wo 



riL20 wo 



1 



Also, 



ii0j20'-^;30«40'-'i20ii0'-^;30«40"-'i20i20 5 



^LiiO ^Li«0 ^L2ro ^-L2«0 

r ^ '-^«l0i20il0'-'«l0«20j20'-"i30«40jl0'-^«30U0i20 

2^ 2^ iiiii2iii-"iu 2^ ^7^^ 

3\j2h---li L1L2 



Zl...Z4 



L1...L4 



^Li«0 i^LalO ^L2l0 /^LilO 
'-'h0l20h0'-'h0l20j20^h0h03i0^h0h0320 ^Lg 



< 



\ hhjiLi 



'-^Zi0«20ji0 I 



2/5 



' X] ^'l^'2ri2 {<^Z^oi20i2o} 
hl2j2Ll 



X] ^'arur^^ {^/touojio} w X] ^'sr^r^j |Cj3ol40i2o} 



< 



\ hkjiLi 
To sum up, we have obtained 

m 



\'^«i0i20ji0/ 



hli32L2 



L 



2/5 



X ^'i^'2ri2|C;iW20j2o} 



3/2 



hh32Li 



< 



{2i+iy 



X (EI 



(2ITTF pS«l,i2,«3 ^'ir'2rZ3 {^^oi20i3o} 



2/5 



< 



Sii;2iiii ^^i^'2rji {C/^ol20jio} /-^i 

By an identical argument we obtain also 



1/2 



< 



^hhjiLi^h^h^ji {^/iW20jio} /jl 



1/2 



(79) 



(80) 



(81) 



(82) 



(83) 



Now we can adopt exactly the same line of reasoning as in the proof of Proposition [9l so that by 
trivial manipulations we deduce that (j7ip and (|72p are indeed sufficient to have that the RHS 
of (I52D and (1551) converges to zero as / — > +00. □ 
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5.4 The case of a general q: results and conjectures 

The following proposition gives a general version of the results proved in Sections 15.21 and 15. 3[ 
The proof (omitted) is rather long, and can be obtained along the lines of those of Proposition 
[9] and Proposition 1101 

Proposition 11 Fix g > 4. Then, a sufficient condition to have the asymptotic relation ^ 
in the case p = q — 1 is the following; 



lim < sup — h sup 



= lim (sup ^^^4^^ + sup 1=0. (84) 

Remarks. (1) As in the proofs of Proposition [U] and Proposition 1 101 a crucial technique in 
proving Proposition [11] consists in the simplification of sums of the type 

'-'hm-^l2m2^ L-^Milsm^'-' L2M2jini'-'hmil2rn2^ LsMshms^ L4,Mij2n2^ 

mim2m:j 
Ml. ..Mi 

by means of the general relation 

Ef^LxMx ^LsMs _ ^LscMi /of.\ 
'-'hmil2m2^limil2m2 ~ °Li°Mi- \°'^) 



mim2 



This basically means that, if in (j85p each Clebsch-Gordan coefficient is represented as the vertex 
of a connected graph, then it is possible to "reduce" such graph by cutting edges corresponding 
to 2-loops - see [11] for a more detailed discussion on these graphical methods. 

(2) Note that, since g > 4 and according to Part C of Theorem [S] condition (I39p is only 
necessary to have the CLT ([^0]) . so that ([84]) cannot be used to deduce the asymptotic Gaus- 
sianity of the frequency components of Hermite-subordinated fields of the type Hq [T]. Some 
conjectures concerning the case q>4:, p^q — 1 are presented at the end of the section. 

(3) Observe that, in terms of the random walk {Zn} defined in (|57p - (|58p . 

'''b.^ ^'''^ = F{Zq_i = X\Zq = l} and 'i'^ ^ = P{Zi = X\Zq = l}. 

^q,l Tg.Z 

As mentioned before, the relation (|39p (which implies (|40p). in the general case where g > 4 
and p q — 1, is still being investigated, as it requires a hard analysis of higher order Clebsch- 
Gordan coefficients by means of graphical techniques (see for instance |28l Ch. 11]). At this 
stage, it is however natural to propose the following conjecture. Recall that we focus on the 

CLT (HO]) because of the equality in law t/''^ (x) = ^j^^aiQ-^q, and Corollary H 
Conjecture A ( Weak) A sufficient condition for the CLT is 

r pM q+l-p,l;X , „^ 

lim max sup (87) 



-F-'J A l^L P,Lq-2'^ q+l-p,l;L 

lim max supPjZp = X \ Zq = 1} = . 

Z— >oo l<p<g— 1 A 
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It is worth emphasizing how condition (j87p is the exact analogous of the necessary and 
sufficient condition ([50]) . estabhshed in [T7] for the high-frequency CLT on the torus T = [0, 2tt). 
This remarkable circumstance may suggest the following (much more general and, for the time 
being, quite imprecise) extension. 

Conjecture B (Strong) Let T be an isotropic Gaussian Eeld defined on the homogeneous 
space of a compact group G, and set T^'^^ = Hq (T) (q > 2). Then, the high-frequency compo- 
nents of T^'^) are asymptotically Gaussian if, and only if, it holds a condition of the type 

1- ^ p,xl^ q+l-p,l;X , . 

lim max sup — — = U , (88) 

where G is the dual of G, Iq is some point at the boundary of G, and the convolutions T and T* 
are defined (analogously to ^U2 \) -[U5 \) ) on the power spectrum of T, by means of the appropriate 
Clebsch-Gordan coefEcients of the group. 

We leave the two Conjectures A and B as open issues for future research. 

Remark. (On ^^no privileged path" conditions) In terms of Z, condition (|87p can be further 
interpreted as follows: for every I, define a "bridge" of length q, by conditioning Z to equal / 
at time q. Then, (j87p is verified if, and only if, the probability that the bridge hits A at time q 
converges to zero, uniformly on A, as Z — > +cxd. It is also evident that, when (j87p is verified for 
every p = 1, <? — 1, one also has that 

lim sup P [Zi = Ai, = Ag_i \Zq = l]=0, (89) 

l^+oo Ai,...,Ag_i6N 

meaning that, asymptotically, the law of Z does not charge any "privileged path" of length 
q leading to /. The interpretation of condition ()89p in terms of bridges can be reinforced by 
putting by convention Zq = 0, so that the probability in ([89|) is that of the particular path 
^ Ai — > ... — > Ag_i I, associated with a random bridge linking and I. 

6 Further physical interpretation of the convolutions and con- 
nection with other random walks on hypergroups 

6.1 Convolutions as mixed states 

We recall that, in quantum mechanics, it is customary to consider two possible initial states for 
a particle, i.e. those provided by the so-called pure states, where the state of a particle is given, 
and those provided by the so-called mixed states, where the state of the particle is given by a 
mixture (in the usual probabilistic sense) over different quantum states. We refer the reader to 
[13j for an introduction to these ideas. From this standpoint, the quantity F^; defined in (j54p 
is the probability associated to a mixed state, where the mixing is performed over all possible 
values of the total angular momentum. To illustrate this point, we use the standard bra-ket 
notation \10) to indicate the state of a particle having total angular momentum equal to / and 
projection on the z-axis. By using this formalism, the quantity Fg ; can be obtained as follows: 

(i) consider a system of q particles ai, such that each aj is in the mixed state S according 
to which a particle is in the state |A;0) with probability Tk/T^, {k >0); 
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(ii) obtain Tg^i as the probability that the elements of this system are coupled pairwise to form 
a particle in the state \10). 

Now denote by ^p^\\o) the event that the first p particles ai, ...,ap have coupled pairwise to 
generate the state |A0). Then, one also has that 



^ P+^M q-p,l;X 



Pr {the q particles generate |/0) | ^p,\xo) }• (90) 



1,1 



In particular, relation (j90p yields a further physical interpretation of the "no privileged path 
condition" discussed in 



6.2 Other convolutions and random walks on group duals 

Random walks on hypergroups, and specifically on group duals, have been actively studied in 
the seventies - see [S] Ch. 6]. Our aim in the sequel is to compare our definitions with those 
provided in this earlier literature, mainly by discussing the alternative physical meanings of the 
associated notion of convolution. We recall from Section [3] that, starting from the Wigner's 
D-matrices representation of S0{3), we obtain the unitary equivalent reducible representations 
{D^^{g)(g)D^^{g)} and {®\''+l^_^^^D\g)}. Now note Xl{g) the character of D\g); for all g G 50(3), 
we have immediately 

h+h 

Xh{9)Xh{9) = X] • 

l=\l2-h\ 

In [SI p. 222], an alternative class of Clebsch-Gordan coefficients {C\^i^]^q : /i, ^2, ^ > 0} is defined 
by means of the identity 



21 

which leads to 



^W^IG " (2/i +^1)12/2 + 1) ^^'^'^^ 



where we use the same notation as in [28J and in many other physical textbooks, i.e. we take 
to represent the indicator function of the event I/2 — ^i| < ' < ^2 + ^i- Of course 

'2+'l 21 -\- I 

l=\l2 — ll\ 

As observed in [8], relation ([9T]) can be used to endow 50(3) with an hypergroup structure, 
via the formal addition h + I2 — ^^/i/jIG" i^i • ^ ^ 0} be a collection of positive 

coefficients such that T/ = 1. The convolutions and *-convolutions of the {Ti} that are 
naturally associated with the above formal addition are given by 

^2,1 = "^Ti^Fi^Cli^^fj , Ts^z = X r2,Lirz3C[^;g|c., ... (92) 

Ll,lq 
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and, for p > 2, 

h Ip Li...Lp_2 

As shown in [8j, the objects appearing in ()92|) -(|M1) can be used to define the law of a random 

walk Z = {Zn : n > 1} on N (regarded as an hypergroup isomorphic to SO (3)), exactly as we did 
in ([57j) - ([58|) . In particular, since = ^ = 1, one has that T* ^.^^ = F{Zp = l\ Zi = li}. Also, 

the convolutions (|92p - (j94p (and therefore the random walk Z) enjoy a physical interpretation 
which is interesting to compare with our previous result. To see this, assume we have two mixed 
states and : in state , the particle has total angular momentum li and its projection on 
the axis z takes values rn-i = — /i, /i with uniform (classical) probability (2/i + 1)^^; analogous 
conditions are imposed for H^j- Let us now compute the probability Pr{/ | H^^jH^j} that the 
system will couple to form a particle with total angular momentum / and arbitrary projection 
on z. Start by observing that the probability that a particle in the state l^imi) will couple 
with another particle in the state \l2m2) to yield the state \lm) is exactly given by {Cf^^i^^^}"^ ■ 
Hence, with straightforward notation, 

= ^ Pr{/ I l^imi) , |/2"22)}Pr{mi,m2} 

mim2 

= EP'-{'ll'.-).l'^-^» 2i.Vl2faVl 

mim2 

l^hnnhm^j 2/1 + 12^2 + 1 



m m\m2 



^ {hhl} 1 _ 2/ + 1 {hhl} _ 
^2/1 + 1 2/2 + 1 2/1 + 12/2 + 1 '^'"''^ ' 



m 



It follows from (|95p that the quantity Fg / can be obtained as follows: 

(i) consider a system of q particles ai, such that each Uj is in the mixed state H according 

to which a particle is in the state \ku), u = —k, ...,k, with probability {2k + l)~^Tj./r^ 
{k > 0); 

(ii) obtain Tq^i as the probability that the elements of this system are coupled pairwise to form 

a particle in the state \lm) , any m = —I, ...,/. 

To sum up, both convolutions F and F can be interpreted in terms of random interact- 
ing quantum particles: F-type convolutions are obtained from particles in mixed states where 
the mixing is performed over pure states of the form \kO); on the other hand, F-type convo- 
lutions are associated with mixed state particles where mixing is over pure states of the type 
{\ku) : u = —k, k}, uniformly in u for every fixed k. 



7 Application: algebraic/exponential dualities 

In this section we discuss explicit conditions on the angular power spectrum {C; : / > 0} of the 
Gaussian field T introduced in Section [21 ensuring that the CLT (HOll may hold. Our results 
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show that, if the power spectrum decreases exponentially, then a high-frequency CLT holds, 
whereas the opposite implication holds if the spectrum decreases as a negative power. This 
duality mirrors analogous conditions previously established in the Abelian case - see [T7j. For 
simplicity, we stick to the case q = 2. Note that the results below allow to deal with the 
asymptotic (high-frequency) behaviour of the Sachs- Wolfe model ([B]). 



7.1 The Exponential case 

Assume 

Q « (/ + l)°exp(-/) , a>0. (96) 

To prove that, in this case, (j40p is verified for g = 2, we will prove that (j66p holds (recall the 
definition of Ti given in (I5ip ). For the denominator of the previous expression we obtain the 
lower bound 

oo [2i/3] 

^ Tij^i^icf^Qi^of > ^ r/ir/_ii(Q%;_,^o)^ 

h,h=l h=[l/-i\ 

[2V3] 

« exp(-/)/2("+i) iCl%_,^,f (97) 

h = [i/3] 



and in view of [28], equation 8.5.2.33, and Stirling's formula 



h=[l/S\ 

I 



[2V3] ,21+1 



exp(-o/2(°+^) y 

Zi = [i/3] 

/(2/i)2'l+V2(2/_2/i)2'-2Zi+l/2\ 



(2Z)2'+l/2 
[2V3] 



;l/2 



« exp(-0/'(°+^) y «exp(-/)/2(°+i)/V2. 

Zi = [V3] (i-'l)'^^ 

On the other hand, recall that by the triangle conditions (Section [3]) {Cf^^^i^,^'^ = unless 
^1 + ^2 Hence 

snvY^Ti.Ti, {QWor <i^supexp(-/)/?+i 

f oo ^ 

/_/,p+i + ^exp(-n) |/i+nr+H «exp(-0/'^"+')/'/^ 



X 

I u=l 

It is then immediate to see that that (I66D is satisfied. 
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7.2 Regularly varying functions 

For q = 2, we show below that the CLT fails for all sequences Ci such that: (a) Ci is quasi 
monotonic, i.e. Q+i < C;(l + K/l), and (b) Ci is such that liminf;__>oo > 0- 

particular, a necessary condition for the CLT (I40p to hold is that Ci/Ci/2 0. This is exactly 
the same necessary condition as was derived by [17J in the Abelian case. For the general case 
q > 2, we expect the CLT fails for all regularly varying angular power spectra, i.e. for all Ci 
such that liminii^oo Ci/Cai > for all a > 0. Note that we are thus covering all polynomial 
forms for Cf^. 

Since (|66p only provides a sufficient condition for the CLT, we need to analyze directly the 
more primitive condition (|39p for m = (however, the case m ^ entails just a more complicated 
notation). We consider first an upper bound for the square root of the denominator of (j39p . 
which is given by . 

We have 



a 



(2) 



(2j i + l)(2j2 + 1) 
47r(2/ + 1) 



c: 



10 



ii0j20 



ii0i20 



Jl J2 



Att{21 + 1) 

1 °° 2 

ii 02= ji 

1 f 1 °° 2 

^ ^E^^-i(2ji + l) sup C,, <i^Q/2 

I^i2>iiji+i2>« J ^2=0 

where we have used the relation ^fj" (Cj lOjjo)^ ~ (^jiozo)^' ^^^^ 



«2 + «l 



sup 

j2>jl,jl+j2>l 



C,, < , and ^ Q\%,o =1 



Z=|«2-«l| 

For the numerator of ()39p one has that it is greater than 

(2ii + l)(2j2 + l) 



E ^h^h' 



Jl J2 



(47r(2/ + l))2 



^ (2ji + l)(2j2 + l) 

- (47r(2/ + l))2 



Jl J2 



Er* /o; _Ln/^w /-y/o 
'-^«il^fi-i-ijL'iiOjio'-^«iOiiO'-"iiOi2o'-'iiOi20 

56*2 I ^20^-^ 0^*20^20} 



> C 



(4vr) 



5C72 <^ a 



10 y 



>KC^ 



The left-hand side of condition ([39]) is then bounded below by lim/^00 (-f^iCf) /(-^2C'f/2) ^ 
so that the CLT (iOD cannot hold. 
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